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Abstract 

After a brief review of previous work, two exactly solvable two-dimensional mod- 
els of a finite Coulomb fluid in a disc are studied. The charge correlation function 
near the boundary circle is computed. When the disc radius is large compared 
to the bulk correlation length, a correlation function of the surface charge density 
can be defined. It is checked, on the solvable models, that this correlation func- 
tion does have the generic long-range behaviour, decaying as the inverse square 
distance, predicted by macroscopic electrostatics. In the case of a two-component 
plasma (Coulomb fluid made of two species of particles of opposite charges), the 
density correlation function on the boundary circle itself is conjectured to have a 
temperature-independent behaviour, decaying as the -4 power of the distance. 

LPT Orsay 02-20 
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1 Introduction 



This paper is dedicated to Jean-Pierre Hansen on the occasion of his 60th birthday. 
Although Jean-Pierre and I have coauthored only one paper [jl|, some 30 years ago, we 
share a longstanding interest for Coulomb fluids. Here is a contribution to this thriving 
domain of research. It consists of a far from exhaustive minireview of previous work on 
exactly solvable two-dimensional models of Coulomb fluids, followed by an original part 
about a finite Coulomb fluid in a disc. 



1.1 A brief review 

The classical (i.e. non-quantum) statistical mechanics of some two-dimensional models 
of Coulomb fluids is exactly solvable. These models have an intrisic interest: this is the 
only case of solvable models for a continuous (i.e. not on a lattice) fluid, in more than one 
dimension. Furthermore, these models can be used as a testbench for a variety of generic 
properties of Coulomb fluids. For mimicking three-dimensional Coulomb fluids in two 
dimensions, one must use the two-dimensional Coulomb interaction which is logarithmic: 
the interaction energy between two unit charges at a distance r of each other is chosen as 
ln(L/r), where L is an arbitrary length which only flxes the zero of energy. Indeed, this 
interaction is the two-dimensional solution of the Poisson equation 

Aln- = -27r5(r). (1.1) 

r 

Of course, these logarithmic models do not describe "real" charged particles, such as 
electrons, conflned in a plane, which nevertheless interact through the three-dimensional 
Coulomb potential 1/r. 

Two models are of special interest: The one-component plasma (OCP), or jellium, is 
made of one species of particles of charge q, embedded in a uniformly charged background 
of the opposite sign. The two-component plasma (TCP), or Coulomb gas, is made of 
two species of particles, with opposite charges ±q. At the inverse temperature /3, the 
dimensionless coupling constant can be chosen as P = jSq'^. It has the remarkable property 
of being independent of the density, and this results into a simple equation of state 
for both models 

f3p = {l-l)n (1.2) 

where p is the pressure and n is the density. 

In the case of the TCP, the model (with pure Coulomb interactions) becomes unstable 



againt the collapse of positive- negative pairs for P > 2, and the equation of state (|1.2| ) is 
valid only for P < 2. If some short-range repulsion (for instance hard cores of diameter 
a) is introduced, the collapse at P = 2 is suppressed and the temperature can be further 
lowered. For small values of the dimensionless density na^, the famous Kosterlitz-Thouless 
phase transition of inflnite order occurs, between a conducting high-temperature 

phase and a dielectric low-temperature phase, at a density-dependent temperature with 
a corresponding P close to 4. 

Finite-size two-dimensional Coulomb fluids exhibit universal properties related to 
conformal invariance. This will not be reviewed here. 
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For the special value F = 2, the OCP and TCP [|T^ are fully exactly solvable: 



The thermodynamic functions and the correlation functions are obtainable. 
For the OCP, it is convenient to define the Ursell function U as 

f/(r,r') =n(2)(r,r') -n^ (1.3) 

where n is the one-body density and n*^^^(r, r') is the two-body density. At F = 2, in the 
canonical ensemble, the OCP maps on a system of free fermions in a magnetic field and 
one finds 

U (r, r') = —71^ exp(— 7rn|r — r'p). (1-4) 

This correlation function has a fast Gaussian decay on a lengthscale of the order of 
(t™)^^/^. Several generic sum rules can be checked. The zeroth and second moments 



of U obey the two Stillinger-Lovett sum rules |Tl| which express perfect screening of a 
charged particle of the fluid and perfect screening of an added infinitesimal test charge, 
respectively. The fourth moment obeys the general expression derived, in terms of the 
compressibility, by Vieillefosse and Hansen |T^, 



n 



/7rF\ 
\4 



)7c/(r,.')k-r'rd.= -/3(|)^ (1.5) 



and the sixth moment obeys a recently derived expression 

'7rF\3 r 3 



— j yf/(r,r')|r-rrdr=-(F-6)(8-3F). (1.6) 



For the TCP, the Ursell functions U depend on the species of the two particles which 
are involved. They must be defined as 

UsAr,T') = n^3{r,T')-nsns' (1.7) 

where s,s' = ±1 denote the signs of the particles, n+ = = n/2 are the one-body 
densities of each species, n the total one-body density, and nfj, (r, r') is the two-body 
density for particles of species s and s', respectively. At F = 2, in the grand canonical 
ensemble, the TCP maps on a system of free fermions. A control parameter is a properly 
rescaled fugacity m, which has the dimension of an inverse length. Although the density 
diverges at F = 2, the Ursell functions remain finite. They are 



' 2 " 

m 



2 
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[/,,(r,r') = - — [K,{m\v-v'\)f (1.8a) 



' 2 ' 

m 
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f/s-s(r,r')= j [Kr{m\v-v'\)Y (1.8b) 

where Kq and Ki are modified Bessel functions, which have an exponential decay; the 
correlation length is l/(2m). The short- distance behaviours of these Ursell functions 



confirm general predictions of Hansen and Viot The pair distribution function for 



particles of opposite signs, dominated at short distances by the Bolzmann factor of the 
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Coulomb potential, should behave like |r — r'|~^; the |r — r'|~^ behaviour of Us-s fits with 
this form. However, although the pair distribution function for particles of the same sign 
should behave like |r — r'|^ for T < 1, this repulsive behaviour has been predicted to be 
weakened into |r — r'p"'" for 1 < F < 2. because of screening by a third particle; the 
logarithmic behaviour of the Bessel function Kq in Ugs can be considered as a limiting 
case when F = 2. 

The Ursell functions for the TCP can be combined into a charge correlation function 

U, = 2q'iU+^-U^^) (1.9) 
and a density correlation function 

[/„ = 2(f/++ + [/_+). (1.10) 

Several generic sum rules can be checked on these correlation functions. The zeroth and 
second moment of the charge correlation function ( |1.9| ) obey the Stillinger-Lovett sum 
rules. The zeroth moment of the density correlation function ( |1.10 ) obeys the usual 



compressibility sum rule, while its second moment obeys a recently dicovered sum rule 
[l^,|T6| which seems to be specific to the two-dimensional point-particle TCP: 

f/„(r,r')|r-r'pdr = . — , (1.11) 

^ ' ^' ' 127r[l - (F/4)2] ^ ^ 

Up to two years ago, there were exact results only at F = 2. Then a major breakthrough 
occurred: Samaj et al. [|1^] succeeded in deriving the thermodynamic properties of the 
TCP in the whole range of stability of the model, F < 2. This was achieved by using a 
mapping on the sine-Gordon field theory and results known for the latter. The surface 



tension along a wall could also be obtained in the cases of an ideal conductor wall [18 



and of an ideal dielectric wall [1^. However, although results have been obtained for the 
asymptotic behaviour of the correlation functions [0,0], there are no simple expressions 
for these correlation functions. In the following, we shall consider only the simple case 
F = 2, and the weak-coupling (high-temperature) limit F — > 0. 

At F = 2, many exact results are available for the OCP and the TCP in a variety 
of geometries with walls or on a curved surface. These results will not be described in 
detail, the interested reader might look on cond-mat for Forrester, Jancovici, Tellez .... 
Here, we only review the case of a Coulomb fiuid living in a half-plane bounded by a 
rectilinear hard wall. This hard wall can be taken as the y axis, and the fiuid is supposed 
to occupy the half-plane x > 0. The position r of a particle is defined by its Cartesian 
coordinates x and y. Now, in the definitions of the Ursell functions, it must be understood 
that the one-body densities are functions of the distance to the wall, i.e. in (|1.3|) n"^ must 
be replaced by n{x)n{x'), and in ( |1 . 7| ) ngHg' must be replaced by ns{x)ns'{x'). The Ursell 
functions depend on x, x', and \y — y'\- 

For this half-plane geometry, a generic behaviour of the charge correlations near the 



wall pT|^3| results from the assumption that the fiuid is a conductor obeying the laws 
of macroscopic electrostatics. The charge correlation function along the wall now is long- 
ranged, with only an algebraic asymptotic decay 

U,{x,x',\y-y'\)-P^, (1.12) 
[y - yr 
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when \y — y'\ is large compared to the microscopic scale (the bulk correlation length). 
f{x,x') is a function which is localized near the wall (it has a fast decay as x or x' 
increase beyond the microscopic scale), and / obeys the sum rule 



/•OO nOO 

(3 dx dx'f{x,x') 
Jo Jo 



1.13) 



( |1.12| ) and ( |1.13| ) can be reexpressed by writing that there is a surface charge density cr(y) 
with a correlation function obeying 



(3 < cT{y)cT{y') >-- 



1 



2n^{y-y'Y 



;i.i4) 



At r = 2 and in the limit F ^ 0, along a hard wall, the charge correlation functions q^U 
of the OCP 1^ and Up of the TCP |0 do have an asymptotic behaviour in agreement 

with (ini) . 

For this half-plane geometry, it has been observed [0 that the density correlation 
function of the TCP near the wall is also long-ranged and that ?7„(a; = 0, x' = 0, |?/ — y'|) 
has the same asymptotic behaviour at F = 2 and as F — 0: 



2'R\y-y'Y' 

It is tempting to conjecture that (|1.15|) is valid at any temperature. 



(1.15) 



1.2 The disc geometry 

We now come to the original part of the present paper. In the above, it was always found 
that the laws of macroscopic electrostatics about charge correlations were satisfied by the 
exactly solvable two-dimensional models, when the lengths under consideration are large 



compared to the microscopic scale. Recently, however, a counterexample was found ||24 
A short-circuited circular capacitor was considered: A two-dimensional Coulomb fluid fills 
a disk of radius Ri and the infinite region outside a concentric circle of larger radius R2-, 
while the annulus between Ri and R2 is empty, and the two filled regions are allowed to 
freely exchange charged particles. This geometry is exactly solvable at F = 2 for both 
the OCP and the TCP. The charge Q on the inner disk fluctuates. Even for macroscopic 
values of -Ri, i?2, with R2/R1 > 1, it was found that the variance < > — < Q >^ differs 
from the value predicted by using linear response theory and macroscopic electrostatics 

/3(<g2>_<g>2)= 1 . (1.16) 

Ri 

This, at first sight surprising, disagreement can be explained, on second thought, when it 
is noted that the fluctuations ( |1.16|) involve only a small number of particles, because, at 
r = 2, {< > — < Q >^)/g^ is of order unity. Although, for one disc alone, i.e. in the 
limit i?2 — ^ C)0, no charge fluctuations are found for the solvable models (except in one 
very special case), in agreement with the limit -R2 ^ 00 in ( |1.16D , the failure of ( p..l(j| ) for 



R2/R1 finite casts a reasonable doubt on the predictions of macroscopic electrostatics in 
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the disc geometry. The subject of the present paper is to investigate the surface charge 
correlations of the OCP and TCP in a disc, at F = 2. It will be shown that macroscopic 
electrostatics does hold for the present problem. 

Thus, we consider a Coulomb fluid in a disc of radius R, bounded by a hard wall. 
It is convenient to put the origin at the centre of the disc, and to use polar coordinates 
{r,(p). The disc may be either insulated, with for simplicity a vanishing total charge, or 
grounded. In both cases, from linear response theory and macroscopic electrostatics 
one finds that the correlation function of the surface charge density is given by 

< a(^).(^') >= (1.17) 

where 6 = ^p — ip' is the angular distance between the two points. Actually, the correlation 
function in ( |1.17| ) is defined from the microscopic charge Ursell function Up{r,r' ,\6\), 
which, for R and R\6\ large compared to the microscopic scale, is localized near the 
boundary circle (it has a fast decay as x = R — r or x' = R — r' increase beyond the 
microscopic scale). The relation is 

/■OO noo 

<a{^)a{^')>= dx dx'Up{r,r',\e\). (1.18) 
JO Jo 

Also, possible oscillations as a function of ^ on a microscopic scale are washed out in the 



definition of the surface charge density, ( p. . 1 7|) is a generalization of the rectilinear wall 



case (|1.14| ), which is retrieved from ( |1.17| ) in the limit i?— >^oo,6'— >0, ata fixed value of 
R6 which becomes y — y' . 

In section 2, it will be shown that ( |1.17D is obeyed in the case of the OCP at P = 2. 
In section 3, it will be shown that (|1.17|) is obeyed in the case of the TCP at P = 2. 
A generalization to the disc geometry of the density correlation function on the wall for 



the TCP ( |1.15| ) will also be presented. The high-temperature limit P — > will be briefly 



reviewed in section 4. 



2 One-component plasma at F = 2 in a disc 

The canonical ensemble is used. There are A^ particles of charge g in a disc of radius 
R. Thus, the average density is n = N/{'kR'^). A uniformly charged background insures 
that the total charge vanishes. It is convenient to choose the unit of length as R/ \/N. In 
these units, nn = 1. The general formalism for the OCP in a disk at P = 2 expresses 
the density Ursell function f/(r, r', \9\) = n^'^^{r, r', \9\) — n{r)n{r') in terms of an auxiliary 
function 

N-l I 



as 



where w = rr' exp{i6) and 7(Z + 1, A^) is an incomplete gamma function p5| 

U{r,r',\e\) = -n^exp{-r^ -r'^)\K{w)\\ (2.2) 



We are interested in the behaviour of ( |2.2| ) when r and r' are close to R and R9 is much 
larger than the microscopic scale ri"^/^, i.e. for A^ large and a non-zero value of 9. 
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A related problem has been previously solved by Clioquard et al. p6[. They considered 
the simpler case when the confining circular background extends well beyond the circular 
blob formed by the particles. Then, in (2J.), the incomplete gamma function is replaced 
by the complete one r(/ + 1) = l\ and the sum itself can be expressed in terms of another 
incomplete gamma function, from which the surface behaviour of U could be obtained 
(the surface behaviour of U has also been obtained ||2^ in the case of an elliptical blob 
of particles). In our present case, since for r and r' close to R the sum in (|2.1|) can be 



seen to be dominated by values of / close to A^, we use the asymptotic expression 
7(/ + l, /) ~ (1/2)Z!. Thus our K is just twice the one of Choquard et al. Although we have 
not been able to make a rigorous mathematical proof, the validity of this procedure has 
been checked on Mathematica, in the cases r = r' = R, 6 = n /2 or 6 = tt , hy numerical 
evaluations of a sum equivalent to (|2.1|) . 



Following Choquard et al., we can now express our A" as 1)25 

^-1 o 



2E 

1=0 



exp(w)r(A^, w) 



{N-l)\ 

where r{N,w) is the incomplete gamma function 

T{N,w) = 

in the present case oi \w — N\ ^ ^/N and A^ 



The asymptotic form of 
obtained by rewriting 



as 



T{N,w 

and integrating by parts, with the result 
exp(w) 



dt- 



N -1-tdt 



it 



(AT-l) 



■T{N,w) 



w 



N 



(2.3) 

(2.4) 

oo, can be 

(2.5) 
(2.6) 



in agreement with an asymptotic expansion by Tricomi (see ||2^). In the case of two 
particles on the circle of radius R, since in our units R"^ = N, w = Nexp{i6), and from 
(gj) and one finds 



NatN 



2e-^^A^ 



4Ne 



Using Stirling's formula for A^! in 



gives 



2 \'/^ 



4Ne 



\nNj 



- 1 



(2.7) 



(2.8) 



Using ( p.8|) in (|2.2| ) and reestablishing an arbitrary unit of length such that vrn is no longer 



1, gives the correlation function on the circle 



U{R, R,9) 



-n- 



(2/?sinf 



(2.9) 
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In the case when the two particles are at small distances x and x' from the boundary, 
minor modifications after ( p.6|) give, in the large- limit 



U{r, r\ 6) ~ —n 



2 exp(— 27rnx^ — 27mx'^) 



(2R 



sm 



(2.10) 



Integrating Up = q^U with respect to x and x' does give a correlation function of the 
surface charge density of the form (|1.17|) (here f3q^ = 2), in agreement with macroscopic 
electrostatics. 

The case of a rectilinear hard wall |21[] is retrieved from ( |2.1CI| ) by taking the limit 
i?^oo,6'— i>0, ata fixed value of R9 which becomes y — y' . Then, in agreement with 

, , , ,,, r, exp(— 27mx^ — 27rnx'^) 
Up{x, x\ \y - y'\) ~ -q^n^^ ^ 2.11 



and ( |1.14|) is satisfied. 

In the derivation of the present results, we have replaced 7(/+l, A^) by (1/2)/! without a 
rigorous mathematical justification. However, the validity of (|2.9| ) and ( p.lO|) is supported 
by several checks. First, the correct result ( p^.llp was retrieved in the limiting case of a 
rectilinear wall. Second, using the exact expression ( p.l| ), Choquard et al. have shown 
without approximations that, in the large- limit. 



-i- / drdr'|r-rf f/(r,rO = - 



(2.12) 



and that the surface contribution to ( |2.12| ) is l/(27r); this is in agreement with what is 
obtained by using for U near the surface the expression ( p.lO|) . 



3 Two-component plasma at F = 2 in a disc 

The particles are confined by a hard wall in a disc of radius R. A grand canonical ensemble 
restricted to neutral configurations is used. Both species of particles, of respective charges 
q and —q, have the same rescaled fugacity m. The general formalism |]10| expresses the 
Ursell functions Uss'{r,r', \6\) in terms of Green functions Gss'{Y,r') as 

UsAr,r', \e\) = -ss'm^\GsAr,r')\^ (3.1) 

Because of the symmetry between positive and negative particles, we only need G++ and 
G |_ which, for r, r' < R, are determined by 



,2 



A)G++(r,r') =m5(r-r') (3.2) 



and 

, exp(i(^) ( d i d \ , 
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In infinite space the solution of (p.2|) is [m/ (27r)]ii'o(m|r — r'|). In a disc, it is appropriate 
to use polar coordinates and to write the solution as an expansion of the form 

m °° 

G++(r, if] r', if') = — V [Ii{mr')Ki{mr) + aili{mr')li{mr)] exp[iZ(v5 - ^p')] {r' < r < R) 

ZTl , — 

t=— oo 

(3.4) 

where the first term in the sum is the expansion of [m/ (27r)]i^o(^|r ~ and the second 
term is a "reflected" contribution due to the wall; Ii and Ki are modified Bessel functions 
and a/ is a coefficient to be determined by the boundary conditions. ( p.3|) gives 

m °° 

G-+{r,(p]r',ip') = — [Ii{mr')Ki+i{mr) 

ZTT , 

t=— oo 

—aili{mr')li^i{mr)] exp[i(/ + l){p — ilip'] {r' < r < R) (3.5) 

For r' < R, r > R, m must be replaced by in (|3.2|), and, as functions of r, G++ depends 

only on 2; = rexp(iy9), G ^ depends only on z = rexp(— iyj), and they must vanish at 

infinity. At r = R, and G |_ must be continuous. These conditions impose that the 

/ > terms in G-^-+{R, (p; r', ip') and the terms / < in G—^{R, ip; r', (p') vanish. Therefore, 
the coefficient ai in (|3.4|) and (p.5|) is 



We are interested in the behaviours of ( p. 41) and ( p.5|) when the disc is much larger that 
the bulk correlation length {mR ':$> 1) , 6 = (p — (p' has a fixed (non-zero) value, and r and 
r' are close to R. 

Let us start with the case when the two points are on the disc. In this limiting case, 
when the Wronskian relation [25| Ii{'mR)Ki^i{'mR)+Ii^i{mR)Ki{mR) = l/{mR) is taken 
into account, and after / has been changed into — /, ( ^.41 ) and ( |3.6| ) give 

1 ^ IijmR) 
2'KRfr{Ii-i{mR) 



G++{R, ^; R, ^') = 77^ E T^T^ exp[-i/(^ - ^')] (3.7) 

1=1 



while, when that same Wronskian relation is used, ( p.5|) and (|3.6| ) give 

00 

G.+ {R, ip; R, <^') = E exp[i(l + l)^ - U^')]- (3.8) 



In the large-mi? limit, using in ( p. 71) the large-argument asymptotic expansions of the 
Bessel functions give, up to order l/R"^, 



Ii-i{mR) 2mR 8{mRy 



Although the sums on / which appear in (|3.8| ), and in ( p.7|) when (|3.9| ) is used, are not 
convergent, they can be given a meaning in the sense of distributions; the obvious recipe 
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is to insert a convergence factor p', with \p\ < 1, in each term, and to take the hmit p 
after the summation. One finds, up to order l/R^, 



G++(i?,<^;i?,<^') 



27rR 1 - e- 



-ie 



1 1 + e 



-10 



2mR 1 - e 



+ 



8{mRy 



-1 + 6e-"' + 3e 
(l-e-i'')2 



-2ie 



+ 



where 6 = (f — (f', and 



G^4R,ip;R,ip') 



(3.10) 



(3.11) 



2nR l-e'f>' 

It might be noted that the expression (|3.11| ) is exact, without any large-mi? expansion. 
Finally, using (|3.10|) or (|3.11|) in (Ol) gives 



U++{R,R,\9\) 



m 



sm 



+ 



1 



47r2 (2i?sin| 



+ 



up to order l/i?^, and 



U^+{R,R,\9\) 



m 



sm 



2 • 



(3.12) 



(3.13) 



More generally, if the two points are at small distances x and x' from the boundary. 



using the large-argument asymptotic expansions of the Bessel functions in (|3.4| ) and (3.5) 
(where the first term in the sum, which comes from the expansion of Kq, gives a short- 
ranged contribution which can be omitted) gives, when only the leading term is kept. 



G++(-R — x,ip; R — x' , ip' 



-•m(x+x') 



-W 



and 



G-+{R — x,ip; R — x', ip') 
The Ursell functions become 

U++{R-x,R-x, \e\) 



2'kR 1 - e-i^ 

Q—m{x+x')+iip Y 



and 



U-+{R-x,R-x\\e\) 



2nR 1 - e'^ ' 

^2^-2m{x+x') 

47r2 (2i?sin |)^ 

^2g— 2m(x+x') 

47r2 (^2i?sin| 



(3.14) 
(3.15) 

(3.16) 
(3.17) 



The charge correlation function is 

Up{R -x,R- x', \e\) = 2q^[U++{R ~x,R- x', \e\) - U^+{R -x,R~ x', \e\)]. (3.18) 

Using (|3.16|) and (|3.17 ) in (|3.18|) and integrating it on x and x' does give a correlation 
function of the surface charge density of the form ( |1.17| ) (here Pq"^ = 2), in agreement 
with macroscopic electrostatics. 
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The density correlation function on the wall is 



Un{R, R, \9\) = 2[U++{R, R, \9\) + f/_+(i?, R, \e\)]. (3.19) 



In ( |3.19|) , the leading contribution to f/++ is cancelled by and it is necessary to use 



the full expansion ( p.l2| ), together with (|3.13| ) for obtaining the leading term of Un as 



Un{R, R, \e\) ^. (3.20) 

27r2 (2i?sin " 



2 



( p.20| ) is the generalization to the case of a circular wall of the rectilinear wall result ( |1.15|) . 



A mathematical justification of some of the heuristic steps used in the above derivations 
is given in the Appendix. 

4 High-temperature limit F ^ in a disc 

For both the OCP and the TCP, the high-temperature limit F — is described by the 
Debye-Hiickel theory. The case of a disc has already been investigated by Choquard et 
al. |^0|]. They have shown that, when x = R — r and x' = R — r' are small, the charge 



correlation function isQ 

u,{ryAe\)^-\ (4.1) 

(2i?sin fj 



where k = (2'7r/3g^n)^/^ is the the inverse Debye length. Integrating (|4.1| ) on x and x' does 



give a correlation function of the surface charge density of the form ( |1 . 1 T] ) , in agreement 
with macroscopic electrostatics. 

In the case of the TCP, the density correlation function [/„ is also of interest. Since the 
Debye-Hiickel theory gives a vanishing result (the contributions from the Debye-Hiickel 
Uss and Us-s cancel each other), it is necessary to go beyond the Debye-Hiickel theory and 



to take into account the next term in the renormalized Mayer expansion [|1^, ^ , which is 
proportional to the square of the Debye-Hiickel Up. One finds 

2k(xH-x') 

f/„(i? - X, i? - x', 1^1) -4. (4.2) 

27r2 (2i?sin|) 



On the wall itself, ( [4.2| ) is the same as in the case ( |3.2(J| ) of F = 2. 

5 Conclusion 

Since, in two dimensions, we could have some doubts about the validity of macroscopic 
electrostatics, which is used for deriving the correlation function (|1 . 1 7| ) of the surface 



^Actually, Choquard et al. have given, instead of (4.1), a more general expression involving the 
arbitrary length L in the logarithmic potential ln(L/r). We have argued that, in two dimensions, 
there should be no fluctuations of the total charge on the disc, even in a grand canonical ensemble, and 
that this condition imposes that the limit L 00 should be taken. Then, the expression of Choquard et 



al. reduces to (4.1) 
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charge density for a Coulomb fluid in a disc, we have checked this correlation function on 
exactly solvable models, the OCP and TCP at F = 2. That correlation function had also 



been obtained in the Debye-Hiickel theory |30], i.e. in the weak-coupling F — limit. 



Thus, these calculations indicate that macroscopic electrostatics is valid for describing the 
correlation function of the surface charge density. Incidentally, this correlation function 
(times (3) can be formally written as a Fourier series 

1 I oo 

P < a(^H^') >= -^-j^j^-^ = _g,,os(;»). (5,1) 

This Fourier series has no / = term, another indication that the total charge has no 
fluctuations. 

In the case of a TCP, the density correlation function on the wall has the same form 
( p.20| ) at F = 2 and in the high-temperature limit F — >^ 0. It is tempting to conjecture that 
this form is valid at any temperature, at least in the range F < 2. A similar universality 
of (|1.15| ) had been conjectured 0] in the case of a rectilinear wall. Proving (or disproving) 



these conjectures is an open problem. 



Appendix 

In section 3, large- argument expansions of Bessel functions have been used, although these 
expansions have been inserted into infinite series involving also large values of the index /. 
Furthermore, these series in / have been made to converge only through the introduction 
of an ad hoc convergence factor. In the present Appendix, some justification is given. For 
simplicity, only one case is considered, the calculation of G-+{R, {p;r' = R — x', ip'). 
From ( |3.5|) , (|3^) , and the Wronskian relation. 



1 °° Ti(mr') 



Using the uniform asymptotic expansions of the Bessel functions [^, appropriate here 
since the arguments mr' and mR are large but the index / may also be large, gives after 
some algebra 

G_+(i?,(^;r',y,') ~ ^Eexp[-mx'(l + -^)V2]exp(i/^) (A.2) 



where 6 = ip — ip'. We are interested in the mR oo limit of ( |A.2| ), for a fixed value of 



mx' and 6^0. For obtaining this limit, the Borel summation method is used. The sum 



in (|A.2|) can be written as 

S= I dtexp(-t)/(t) (A.3) 



oo 







where 

oo ;2 4.1 

/(t) =Eexp[-mx'(l + -— )V2]exp(i/^)-. (A.4) 
1=0 
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Since the sum ([A .41) is absolutely convergent, the limit and the sum can be interchanged, 
giving 

. lim /(t) = exp(— mx') exp(te'^) (A. 5) 

mR^oo 



Using the limit ( |A.5| ) in (|A.3|) gives the large-mi? behaviour 

„—mx'+iip 1 

G,4R, ^■,R- a:', V^O - (A.6) 
in agreement with (|3.15|) . 
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